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SUMMARY 

The linearized attenuation theory of NACA 
Technical Note 3375 is modified in the following 
manner : (a) an unsteady compressible local skin- 
friction coefficient is employed rather than the 
equivalent steady-flow incompressible coefficient; ( b ) 
a nonlinear approach is used to permit application 
of the theory to large attenuations; and (c) transition 
effects are considered , Curves are presented for 
predicting attenuation for shock pressure ratios up 
to 20 and a range of shock-tube Reynolds numbers. 
Comparison of theory and experimental data for 
shock-wave strengths between 1.5 and 10 over a unde 
range of Reynolds numbers shows good agreement 
with the nonlinear theory evaluated for a transition 
Reynolds number of 2.5 X 10 s . 

INTRODUCTION 

The increasingly widespread use of the shock 
tube as an aerodynamic testing facility has led to 
the closer investigation of the flows present in 
such tubes. In particular, since the deviation of 
these flows from those predicted by perfect fluid 
theory is often of large magnitude, these devia- 
tions have been investigated fairly thoroughly. 
Several such studies, either of an experimental or 
theoretical nature, may be found in references 1 
to 11. Investigations of the boundary layers in 
shock tubes have been made in some of the afore- 
mentioned references as well as in references 12 
to 17. This list of references does not cover the 
complete field of literature existing on these topics 


hut is representative of the various general 
treatments. 

Consideration of the entire flow field from the 
leading edge of the expansion wave to the shock 
wave is necessary to obtain an accurate picture of 
the waves traveling along the shock tube. These 
waves arc responsible for the deviations from 
perfect fluid flow in shock-wave strength (attenu- 
ation) with distance, in pressure and density at a 
given distance with time, and so forth, which have 
been noted by various investigators. The analysis 
of reference 1 was the first to treat this complete 
flow. Figure 1, a reproduction with minor changes 
of figure 1 of reference 1, is the basic wave diagram 
of the unperturbed shock-tube flow si lowing the 
various flow regions to be considered with typical 
characteristics and particle paths. This linearized 
analysis (red. 1) was based on an averaged one- 
dimensional nonsteady flow in which wall-shear 



Future I. — Distance-time plot of shock-tube flow. 


1 Supersedes NACA Technical Note 4347 by Robert L. Trimpt and Nathaniel B. Cohen, 1958. 
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and heat-transfer effects generated pressure waves 
to perturb the perfect fluid flow. This averaging 
process essentially implied thick boundary layers. 
The expansion wave was treated as a “negative 
shock” or zero-thickness wave. The resulting 
perturbation equations then lunged on the evalu- 
ation of the local skin-friction coefficient which 
in reference 1 was assumed to be given by an 
equivalent incompressible steady flow. Conse- 
quently, the application of the results of reference 1 
was limited to shock pressure ratios in which this 
assumption for Cf was valid, although the analysis 
was still applicable for other pressure ratios when 
the proper choice of c f was employed. The as- 
sumption of incompressibility should apparently 
eliminate the strong shock pressure ratios from the 
range of validity. 

Solutions to the laminar boundary-layer equa- 
tions employing a linear viscosity-temperature 
relation (refs. 2, 3, and 14) show that the nonsteady 
character of the flow is such that the equivalent 
laminar steady-flow assumption is in error, irre- 
spective of compressibility', for most conditions 
except that existing in the cold-gas region a for 
strong shock waves. On the other hand, the 
turbulent boundary layer is not nearly so sensitive 
to the unsteady" character of the flow. Ref- 
erence 15, which assumed a one-seventh-power 
velocity' profile similar to that of reference 1, 
reported that even for infinite shock pressure 
ratios the effect of unsteadiness would produce 
only' a maximum variation in turbulent skin 
friction of 5 percent in the cold gas and of 
22 percent in the hot gas. 

The only- other attenuation analysis to date 
that considers the entire flow field is that of 
reference 2. This analysis is similar to that of 
reference 1 in that it is a small-perturbation 
approach using traveling waves and a negative 
shock, tlie major difference being that the pressure- 
wave generations arise because the boundary- 
layer-displaeement thickness changes with time. 
(The boundary ' -layer-displacement thicknesses of 
ref. 15 are used.) Flows with thin boundary 
layers having a linear viscosit y- tempo rat lire vari- 
ation are required for this treatment to apply. 
The attenuations predicted by references 1 and 2 
for turbulent boundary layers agreed within 
10 percent for shock-pressure ratios up to 6 in 
spite of the marked differences assumed in the 
mechanism for handling the wall effects. The 


perturbations in the flow behind the shock show 
a larger difference between the two approaches. 

The deviations from ideal theory discussed arise 
for the most part from wall effects, that is, the 
perturbations in the shock-tube flow caused hy r 
wall shear and beat transfer. Much recent work 
has been done using the shock tube as a testing 
medium to provide very high- temperature flows 
of short duration. (Sec, for example, refs. 9 and 
18.) In these cases, deviations from ideal fluid 
flow will also arise because the air at high temper- 
atures does not behave as an ideal fluid. It. 
would be difficult to separate the real-gas effects 
from the wall effects; therefore, the present 
analysis, like those of references 1 and 2, is con- 
cerned only with the effects of wall boundary 
lay^er upon the in viscid outer flow, the fluid being 
considered as an ideal gas. 

The turbulent theory" of reference 1 has been 
compared with experimental data for attenuation 
in references 1, 7, 8, and 10 and good agreement 
has been found in general. Predicted pressure 
perturbations in the hot gas by the method of 
reference 1 agreed well with the experimental 
results reported in the same paper. Fair to good 
agreement between theory and experiment is 
reported in references 7 and 10 for the hot-gas 
average density' variation with time in the flow 
behind the shock wave; poor agreement is reported 
for the cold-gas flow where the finite expansion fan 
has been treated as a negative shock. 

Since the deviations from the inviseid fluid flow 
often become large in cases of aerodynamic shock- 
tube testing, the linear, or small-perturbation, 
theories of references 1 and 2 are no longer 
applicable and recourse must be made to some 
sort of nonlinear approach. 

In order to obtain an exact theory for predicting 
the perturbations in a shock-tube flow, a rigorous 
treatment would be required first to the solution 
of the boundary-layer flows. The bound ary-laycT 
equations would have to be solved not only in 
region ($ but also inside and after the expansion fan 
which is considered to be of finite extent. For 
laminar flows the main difficulty would probably 
be the correct handling of the viscosity variation 
across the boundary layer. For turbulent flow a 
rigorous treatment appears to be impossible with- 
out a tremendous increase in knowledge of the 
mechanics of turbulence. Once the boundiiry- 
layer solutions were determined, the vertical 
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velocity at the edge of the boundary layer could, 
if the boundary layers were thin, be used in the 
manner of reference 2 to determine the local 
pressure waves generated. 

The second major difficulty in obtaining a 
rigorous perturbation solution would arise from 
the treatment of the entropy discontinuity. The 
theoretical contact surface increases in extent with 
distance progressed down the shock tube due to 
mixing and diffusion (the former is the major 
influence) . This process not only generates pres- 
sure waves but also alters the reflected and trans- 
mitted wave strengths of the pressure waves 
generated by the boundary layers. 

If a rigorous solution such as that just described 
was available, then it would without question be 
the one to be employed. The theory of reference 
1 and the theory of the present report assume 
that the wall effects can be averaged across the 
flow. Tli is assumption introduces errors because 
no such physical mechanism exists for the instan- 
taneous transmission of these effects across the 
flow. In the absence of the rigorous solution there 
is no evidence to indicate that the errors intro- 
duced b} T a shock-t ube perturbation theory based 


on an averaging process are of a larger magnitude 
than those introduced by the neglect of the 
aforementioned considerations required for a rig- 
orous solution. In addition, there is the possibil- 
ity that an averaging process might be more 
applicable as the boundary layer fills a greater 
part of the shock-tube cross-sectional area. Con- 
sequently, the extension of the method of reference 
1 in (lie present analysis is justified. 

In the present paper the analysis of reference 1 
is first modified to eliminate the restrictions 
imposed by the incompressible equivalent steady- 
flow assumption for local skin friction; and then 
a nonlinear theory is derived which permits appli- 
cation of the analysis to large attenuations. It 
will be assumed that the reader is familiar with 
the basic theory and assumptions of reference 1 
so that repetition in this paper may be avoided. 
This modified theory will be compared with 
experimental data covering a wide range of flow 
variables. The theoretical and experimental stud- 
ies reported herein were conducted at the Gas 
Dynamics Branch of the Langley Laboratory 
during 1955 and 1956. 
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SYMBOLS 

velocity of sound 


constant defined by equation (8) 
local skin-friction coefficient, 2 rJpU 2 
coefficient of specific heat at constant volume 
coefficient of specific heat at constant pressure 

hydraulic diameter, - 

rerun e ter 


functions defined by equation (89b) 


constant defined by equation (36) 


ratio of contributions of P waves to total waves generated in region f3, 
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linear attenuation with first subscript describing boundary layer appro- 
priate to region a and second subscript to region /3; that is, 
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turbulent with n — 7 and region j3 as tran- 
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fixed distance along shock-tube axis 

distance the shock moves from a given point until the effects of transi- 
tion in the flow generated at that point first influences the shock 


reciprocal of velocity exponent In boundary layer 


v_ 

’U 



i_ 

n 


arbitrarily denotes subdivisions I, II, ITT, etc. of hot-gas region 
nonlinear at tent nation with first subscript describing boundary layer 
appropriate to region a and second subscript to region 
Prandtl number 


character istic parameter, 


2 c, 
R 


a± U 


effective characteristic wave parameter defined by equations (43) 
static pressure 

gas constant ; Reynolds number 

. . m* 

Reynolds number of transition, — — 


entropy 

temperature 

wall temperature 

adiabatic wall temperature 

time 

f r e e-s tr earn v elo ci t y 


shock velocity 

velocity of wave which generates flow 
velocity in boundary layer 

distance along shock tube from diaphragm station 


distance from surface 

function defined by equation (1312) 

function defined in equations (31) 

ratio of specific heats, c p /c v ; assumed as 1.40 for computations 

length of segment into which shock tube is divided for nonlinear treatment 


characteristic derivative in boundary layer, 


&(_) ,J s* a( ) 

dx d df 


. . . . . .. . , 3 3 5 I 6* 1 5( ) 

characteristic derivative m boundary layer, ^+^7 -- -p ^ 

boundary-layer thickness; also indicates differential quantity 

boundary -layer displacement thickness, f(-p) dy 

characteristic derivative in potential flow, — + (U±a) 


contribution to attentuation due to P waves in region (3 
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7)r recovery factor, assumed equal to 0.85 for laminar flow and 0.90 for 

turbulent flow 

d boundary-layer momentum thickness, J ™ <ly 

d 0 boundary-layer momentum thickness at £— £ 0 

for £*>£ rf and ifj*^£ rf , its value is 1 

n coefficient of viscosity 

v coefficient of kinematic viscosity 

£ distance flow has progressed along surface 

£ rf distance flow has progressed along surface at entropy discontinuity 

£ 0 distance flow has progressed along surface when 6=6 0 

£* distance flow has progressed along surface at transition 

p density 


<!>P,a } <t>P,&, 4>Q,0 } <Ps,a, 
<f>P 7 4>Q 

or U n 


wall shearing stress 

influence coefficients, defined in appendix A 
influence coefficients defined by equations (44) 
compressibility correction 
exponent in viscosity-temperature law, p ~ T w 

subdivisions of hot-gas region 0 for nonlinear treatment (see fig. 8) 


Subscripts: 

Letter subscripts not included in the symbols 
defined above refer, in general, to values at points 
or within regions shown in figure L Numbered 
subscripts refer to points in figures 8 and 9. 
Exceptions to be noted, however, are as follows: 
m,n refers to velocity profile parameter m,n 
o perfect-fluid value 

t at time t 

rs evaluated immediat ely behind shock, that 

is, point v located at x— U s t 
x at distance x 

T denotes attenuation with transition 

X arbitrary condition in shock-tube free 

stream 

std denotes XACA standard atmospheric 

conditions 

A prime on a symbol indicates a quantity 
evaluated at reference temperature. 

THEORY 

DERIVATION OF EXPRESSIONS FOR LOCAL SKIN-FRICTION 
COEFFICIENT 

The skin-friction coefficient for the flow behind 
wave-induced flows will be found by an integral 
method. An incompressible skin-friction coeffi- 
cient will first be determined and then a simple 
compressibility correction will be applied. 


Tlie integrated equation of motion for the 
incompressible boundary layer with zero pressure 
gradient is (see ref. 14) : 



The form parameter 5*/6 is assumed to be con- 
stant; this has boon shown to be true for the 
unsteady wave-induced laminar flow (sec ref. 14) 
but has not been completely established for 
turbulent flows. Equation (1) then becomes 



Since the resulting expression for c f will ulti- 
mately be used in the attenuation formulas wherein 
the integral J*C/(£M£ is desired, the variable £ is 
introduced. The variable £ is defined as the dis- 
tance a particle in the free stream has moved to 
reach the point ( x,t ) since acceleration by the 
passing wave which originated at r=0 at £=0 and 
which travels with velocity U w . Thus, 

£0, C U w t-x ) 


(3) 
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In tlic case of flow in a shock tube, the value of 
TJw is U s for the flows induced by shock waves. 
If the assumption of reference 1 is followed and 
the expansion wave replaced by a wave of zero 
thickness, moving with the speed of the leading 
edge of the original expansion wave, then 

T <w = c t . 

The differential equation (2) is transformed 
from the x,t coordinate system to the £,£ system by 
using the following derivatives: 



Thus equation (2) becomes 

/M\ /d6\ = JL_ 

\WiU 0 2 E 


eh ere 


..Ii w (l — \\ 


( 5 ) 

CO 


Now equations (2) and (5) arc differential 
equations capable of solution by application of 
the method of characteristics. The slopes of the 
characteristics are 


steady flows, then 




m\ »+ i 


( 8 ) 


For turbulent flows n is the reciprocal of the 
exponent in the fractional power expression 

used to describe the turbulent boundary- 

layer velocity profiles. For laminar flows the 
value of n is one. The C n terms are arbitrary 
constants to be evaluated later and may be com- 
pletely different for the steady and unsteady flows. 

Combining equations (8) and (7) and inte- 
grating yields: 



Ad 1 r (V 
A£ 2E u \p 

0' ,+vA 

(9) 

re JL- c i 

;£rj> 

(10) 


2 - I //+1 


T «-«r 

(11) 

L 1 n+ 1 2E \U 



Substituting equation (11) into equation (8) to 
obtain c f as a function of £ yields: 


[: 


71+1 

C "+ 3 


»+: r " c „ W W 


71+3 

+1 +? ft-fc) 


2 

w+3 


(12a) 


At 1 5* 

Ax U 6 


for equation (2) and 

At 1 5*1 
as U e E 

for equation (5). 

Thus, if the symbol A//Aij is used to denote the 
derivative along a characteristic of slope At/A £, 
equation (5) may be written as 


At? J_ 

A 


( 7 ) 


For steady-flow boundary layers with zero pres- 
sure gradient, it has been established that 
Cf=Cf(&,U,v). If this relation is assumed to hold 
for unsteady flows in the same form as for 


C /({)=^+ 3 



n+3 


+ 


u 



(12b) 


For the special case of 6 0 ~ 0 at £ 0 which cor- 
responds to flow initiation at £ 0 there results 

C/ (t) = F,^» [-({-{.)] " +3 (12c) 

The values of the various F n terms, which arc 
directly related to the hitherto arbitrary C n 
values, will now be determined to match known 
solutions for certain limiting cases. If the ratio 
Uw/U becomes infinitely large, the solution must 
be the same as that for an infinite flat plate in 
contact with a fluid impulsively started from rest 
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at time <=0. Under these conditions 

i-a 0 =ut 

S* 


E=- 


e 


, * /5*\n+3/h 2 /\ n+3 

= ) J 


(•3) 


The solution, known as the Rayleigh solution, to 
sucli an impulsive laminar flow over a plate is 


«/(£) 


= l.,2sf^V' /2 


Q 


\ v J 

2,409 


Rayleigh 


(14) 


Consequently, in order (o match the Rayleigh 
solution for the laminar case (n=l) 




or 


7^ = 0.718 


(15) 


On tlie other hand, if U w jU becomes infinitely 
small, the solution must be the same as that over 
a semi-infinite flat plate in steady flow (that is, 
the so-called Blasius problem). For these con- 
ditions 

757=1 .0 


d«)=c(u^)"' ,+3 j 


(16) 


whereas the corresponding form of the Blasius 
solution for laminar flow is 


t il "1-i/a 

fc)J 

(?) 

\ V / Blasius 


= 2.605 


(17) 


Thus, if the Blasius solution were to be matched, 
it would be necessary for 

7^ = 0.664 

Two possible solutions are then available for 
the laminar incompressible-flow case depending 
upon which limiting value is matched: 

5W>76S 61 2 


Rayleigh limit 


7 7 


( c /)n-l(£) — 0.71s(l 

77 \ 1/2 r~TT “I -1/2 

+ , - 4 fi h-i-e) 0Sa) 

Blasius limit >0^: 

(c/) B -i(f) = 0.604(l 

77 \ 1/2 rrj “ 1 - 1/2 

+ ,005 fv=c) b (i_fe, J (lsb) 


In figure 2 the values for 


v? 


(£—$«>) as de- 


termined by four different means are plotted 
against pressure ratio across a shock in an air-air 
shock tube. The upper branch of the curves 
applies to region (3 behind the shock, and the 
lower branch applies to region a behind the zero- 
thickness expansion wave associated with the 
shock of strength p^ 0 /p m . In addition to the 
values determined from equation (18), there are 
shown values which would be obtained, if the fluid 
were assumed to be incompressible, by the integral 
method of reference 14 and the numerical solution 
to the Prandtl boundary-layer equations. (The 
results of ref. 15 are applicable to this numerical 
solution.) The agreement between both the 
curves of equation (18) and the referenced curves 
is very good. However, since t lie curve based on 
the Rayleigh limit gives a better approximation in 
region 0, which will be shown to dominate the 
attenuation equations, the Rayleigh values of 
7^ = 0.718 and 6*/0=2.4G9 shall be used for t lie 
remainder of this paper. Reference 17 also em- 
ployed a normalized Rayleigh velocity distribution 
in the treatment of flow induced by shock waves. 

For the turbulent case two analogous limiting 
processes are not available in order to determine 
the values of F n . The turbulent boundary-layer 
theory is semiempirical and relies on experiments 
to supply constants for the resulting equations. 
Since no “Rayleigh -type” experiments have been 

U w 

performed, there is no limiting process -jj — 

to apply for the turbulent case. There is, 
however, the semi-infinite flat-plate solution 
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Figure 2. — Plot of incompressible local skin-friction coefficient for wave-induced flows having ratios of wave velocity to 
flow velocity corresponding to those in a shock tube of given shock pressure ratio. The expansion wave has been 
approximated by a negative shock. 1.4. 


corresponding to the limit This solution 


U} 

u 


(sec refs. 1 and 19) assumes a velocity profile 


and results in a skin-friction coefficient 


u 

U \5 

expressed as 


c /(f) — 0.05S1 


[M]~‘ 


(19) 


The combinat ion of equations (16) and (19) results 
in a value of ^7= 0.0581 . 

U /w\ U ^ 

For profiles of the family ^ ^ may be 

shown that 

n 

d~'(n+\)tn+2) 


1 


<5 71+1 

<5* n+2 


6 


n 


( 20 ) 


Therefore, the expression for the onc-seventb- 
powor turbulent velocity profile skin-friction 
coefficient in the nonsteady incompressible flow 
becomes 


<vtt) = 0.0581 (1 +| [7^-u]' 4 


( 21 ) 


Since there was only a minor difference between 
the unsteady-flow values of c f for laminar flow 


" Uw 0 and 


when based on the limiting cases of y 

jj 

^ > 00 9 it is expected that the turbulent-flow 
agreement would be just as favorable if results 


were available for 


lh 


U 


Consequently, equa- 


tion (21) is assumed to be a fairly close approxi- 
mation to the correct answer. Equation (65) of 
reference 15 is very similar to equation (21) but 
was derived in a different way. 
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The skin-friction-eoefficient relation of equation 
(19) corresponding to the on e-seventh -power pro- 
file law is no longer valid at arbitrary large values 
of Reynolds number in incompressible steady flow. 
Instead a logarithmic law is often used. (See ref. 
19.) However, since a power profile is easily 
handled by these methods, the skin friction on a 
semi-infinite plate for these large Reynolds num- 
bers is found to be closely approximated by the 
relation 

C/ =0.0186^^t^-^J _i (22) 

which is compatible with the relation 

If for consistency it is further assumed that 5*/6 
has the value 15/13 (the value for n — 13), then 
the unsteady skin-friction coefficient would be 
given by 

c/i)=001S6(t + l^)*pfcW]-‘ (23) 

No claim is advanced that a 1/13-power profile 
actually exists at these higher Reynolds numbers; 
it is only necessary for equation (22) to be valid 
for steady flow and the value of 5*/d to be 15/13. 

It was shown in reference 14 that the skin 
friction in region a is correctly given by the 
method of characteristics using the boundary 
condition 0—0 for £— 0 at x— — a ( t (that is, on the 
zero-thickness expansion wave) for that part of a 
where x<^Utd/8* and the boundary condition 0=0 
for £— 0 at x= U s t (that is, on the shock wave) for 
that part of a where Ut0/8*gx^Ut. Sketch 1 
shows a boundary-layer momentum character- 
istic. In order to determine the unsteady fric- 
tion coefficient at point c in the region in question, 
UtO/8* ^x ^ Ut, the boundary condition of 0=0 on 
the shock wave (point a) is correct. In the analysis 
which follows, the boundary condition 0=0 on the 
expansion wave (point b) is used instead. Fric- 
tion coefficients are shown qualitatively in sketch 
2. The solid curve represents the case with the 
boundary condition on the shock whereas the 
dashed curve represents the case with the bound- 
ary condition on the expansion wave. In the 
region in question it is seen that the differences 
are not serious. Because of the relatively small 
contribution of region a to attenuation as com- 
pared with the contribution of region £ and also 
in the interest of simplicity, this error is neglected. 



X 

Sketch 1. 



x 

t 


Sketch 2. 
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Tn order U) handle the transitional flows which 
occur behind the waves in a shock tube, some ap- 
proximation for c f in the transition region must be 
employed. Any of a number of assumptions is 
possible in this region. However, in view of the 
many assumptions already present in the theory, 
the least complicating supposition will be em- 
ployed in this case; namely, an instantaneous 
transition is assumed and the value of c s in the 
turbulent region will be taken as the value which 
would be present had turbulent flow existed since 
the initiation of flow. Tn other words, at the 
transition point the local c s changes discon t bili- 
ously from the laminar to the turbulent value and 
the value 0 then applies in both laminar and 
turbulent regions. This assumption was used in 
the logarithmic transitional curve for steady flow 
of reference 19. Figure 3 compares the inte- 
grated skin-friction coefficients of reference 19 
with the curves obtained by the various power 
laws and the foregoing transitional assumption. 
The agreement appears to be very good. 

A simple compressibility correction will be 
based on the intermediate temperature or T f semi- 
empirical method. This correction, expounded 


in reference 20 for laminar boundary layers and in 
reference 21 for turbulent boundary layers, as- 
sumes that the incompressible skin-friction rela- 
tions apply to compressible flow if the properties 
of the compressible flow are taken as some inter- 
mediate value between the wall and free-stream 
values. Thus, if the relation 



applies to a steady incompressible flow, then 



will apply to a steady compressible flow for a cer- 
tain choice of the primed state. The following 
values of the intermediate temperature T' are 
given (see refs. 20 and 21); 

For laminar flow: 

U = i + n.032M 2 +0.58 ( y ~ ] ) ( 26)| ) 
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For turbulent flow: 

^=l+0.035M 2 +0.4o(^-l) (2Gb) 

It may be shown that for an arbitrary reference 
state, wlii eli may be either wall or free stream, 
equation (25) may be put in the form 

/ T\ 7 -^ / '\ — 

" + 3 ( f ) /?+ "( m / (27a) 

c,= F. (~ (27b) 

If a temperature-viscosity relation of the form 



is assumed to apply, then the steady-flow com- 
pressibility correction £l n becomes 






p! 72 


It is assumed that this T f method is also appli- 
cable as a compressibility correction for the un- 
steady-flow skin-friction coefficients. Compressi- 
bility corrections in hot and cold gas are plotted 
against shock pressure ratio in figure 4 for an 
air-air shock tube. The value «=0.8 has been 
used to compute these curves. 

The results of tins section may be summarized 
by the following expressions for skin-friction 
coefficient: 


(<?/) . - 1 =0.7 1 8 (l + 1 .469 f I "'/ •)’((, f ) Qi 

(c/)»- 7=0.0581 (l +f () \ 

w.., ,=0.0186 (. +nr^r)*(7 0 


r (29) 


where 


12 


n 



(28) 




Shock pressure ratio, P^ 0 /P m 


Ficjure 4. — -Compressibility corrections as n function of 
shock pressure ratio for shock-tube flow. 7=1.4; 

w — 0 . 8 . 
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EVALUATION OF THE LINEARIZED ATTENUATION EXPRESSIONS 

The basic linearized attenuation expressions derived in reference 1 are summarized in appendix A. 
The following expressions result from the first identity and subsequent substitution and manipulation 
of the pertinent equations from the appendix: 

Py$ jj8g __ ^gg Pf | P f P c | P c P b | Pb~ P do 

At (If (7* 77 1 (If 


Pf—Pff o _ 1 C id u <rM,—Mf lu 

~ = DJ 0 


oMs—Mg PpQ—Pao ■ 7 1 dap Ppo—Pg 

crM s — P&uPQ&o y ( a Po\ 2 _ i a * 


+ />Jo r /v 


Mg \ 7} P Qpo T 1_ 77 ao J g 0 I ao , 

Ma+rj+l Pao+Qto y a £ 

\77ao/ 

h ^r h -'r>£‘ t ^+d£‘ 


4>S,a 7 




(see ref. 1) the following relation results 


Prs P&o y 77 e Pys J po 

pPo 2 a 0o a ( 


Pis — P d, 

Ppo 


° (2 a~ T)} — r^, CfJ, «7^T" r a ,m "f" Fa, n * C f' a (^3) 


in which the subscript m applies in the laminar range 0 ^£ atj s^£ a , 3 * and the subscript n applies in the 
turbulent range £ a ,p* ^Zaj^Zd- The value of k is defined as £ d /£* for and as 1 for 

The total linearized attenuation is thus made up of the sum of the effects of regions a and 0. 
Th is relation may be expressed mathematically as 


where 


Pvs—PZo f Pv —P»<\ | / Prs—P3<> \ 

Pto V P&o / « \ P&o ) & 

( ! ‘—£r ).“i £»('"■■/. J.t.. c '- j £ 


Now, with C/ expressed in the form 


uo=PnP\r + \^ 


2 _ /T T v—j-/ 

+3/C *.\ «+ 3/ / \ n+3 


2 


integration of equations (35) results in the following equations. (The subscript x designates either 
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region a or region /3.) 


/ P,.—P<h \ ] J r ™ + 3 r 

V M ) x . m .n~2a, 0 D I m+1 


"2 a^T) I m+1 


.(t) m+3 " +I [fe) : «-( 4 r 

xmrwiffmr 


+T„G n a, 




£*Y+ 3 »+ 3 


Rearrangement of terms of equation (37) produces 


/ Pr S —Pn„ \ [ifaJX 1 = y a^pfa 

\ P«, Jx.m,n\jl\V a> J J ZClfoP", 


r m Gja. 


m+3 (V v n \ ™+3/ £ 


mi) 


m + 1 . . m+1 


U \ i 


Tim 




F »+! , v 2 — 1—1 m + 1 p- — 1 n >_ P »+l~I 

(*eL. ter ca ] - °+r te) [• <t r ] 


wlicrc 


'- Jxm \a. vj 


/ g < Pfo m +3 ^ AA m+; * 

Jxra 2ap 0 p m m+1 xmt7xm V/ J 

Note that is a function only of the shock pressure ratio paolpoa and the value of m for any region 
X- The term g xm requires in addition a temperature-viscosity relation which may be either an expo- 
nential type or some other form, such as the Sutherland equation. 

For the case of no transition when the flow is either completely laminar or completely turbulent, 
equation (39a) reduces to 


m+i 2 

v—Tth _ ( l Y ‘+ V a „ p\ m+s 
Po, \I>) \ V m J Jxm x 


7?+l 2_ 

p„s—P(jo/ 1 Y+ 3 /fl'«,7A *+ 3 


o { l \h+3/ a„B\ w+3 

~~\o) \~m) 9xn xn 


The attenuation in this case for a given initial value of ptolp* is obviously dependent only on the two 
parameters expressed as hydraulic diameters of shock-wave travel IjD and shock-tube Reynolds 
number a^Djv^. 

For the case with transition three parameters are required to describe the phenomenon at a given 
initial value of the third parameter is the transition Reynolds number 7? * and enters into 
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the analysis in the following manner: 


^co V v l 
V l ^ £(l 


( 41 ) 


When equation (41) is substituted into equation (39a), the linear transitional attenuation relation 
becomes: 



where the term ^ ~ is a function of Thus, at a given value of p? 0 ly« the linearized atten- 

LJ £ rf 

nation is a function on^ of IJD } and 7?*. 

The attenuation functions g xm , the compressibility corrections Sl xm and their products are presented 
in table I for shock pressure ratios from 1.0 to 20.0 and for m equal to 1, 7, and 13. It may be seen by 
inspection of equation (40) that, forgiven shock tub e 7/7), Reynolds number and no transition, 

the attenuation contribution of the region x is directly proportional to the product (Sl(j) xm . In order to 
demonstrate graphically the behavior of these compressible attenuation functions, they have been 
plotted in figures 5(a) and 5(b). The magnitude of the contribution of region 0 to attenuation increases 
monatomieally with increasing pressure ratio, is always negative in sign, and thus tends to increase 
attenuation. On the other band, for low pressure ratios the contribution of region a tends to increase 



(a) Region a. 



Shock pressure rot io, P^ 0 /P 9 


(b) Region p. 


FinuRE 5. -Variation of compressible linear attenuation functions with shock pressure ratio. y= 1.4; o> = 0.S. 
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attenuation ; above a shook pressure ratio of about 
5.9, the trend is reversed and region a contributes 
compression waves that tend to decrease attenua- 
tion. This reversal is discussed more fully in 
reference 1. 

When no transition is considered, the attenua- 
tion function for the entire flow field is found by 
adding the contributions of regions a and p. For 
the cases where the profile exponent rn is the same 
in both regions, the total compressible attenua- 
tion function has been computed and is shown in 
table I and in figure 6 for values of m of 1, 7, and 
13. This function is given as 

The results of using the methods of references 1 
and 2 are also shown in figures 6(a) and 6(b) for 
values of m of 1 and 7, respectively. 

For the laminar ease the curve in figure 6(a) 
representing the method of reference 1 falls far 
below that of the present report, primarily be- 
cause of the importance of the neglected un- 
steadiness effects as discussed in the introduction. 
The results of reference 2 are also below that of 
the present report (approximately 25 percent for 


i 1 

CM 



(a) Laminar flow; m— 1. 

Figure G. — Plot of total compressible linear attenuation 
functions against shock pressure ratio. 7 = 1.4; w =0.8. 


5G97G3 — G1 3 



(b) Turbulent flow; m — 7. 
Figure 6 . — Continued. 


ro 



(c) Turbulent flow; m= 13. 
Figure 6 . — Concluded. 
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shock pressure ratios from 4 to 10) and show 
better agreement at higher pressure ratios (only 
10 percent below for a shock pressure ratio of 20). 

Agreement between the methods of references 
1 and 2 and the present method is better for the 
turbulent case (m = 7; fig. 6(b)). The neglect of 
unsteadiness has a smaller effect upon the results 
of reference 1, although the effects of compressi- 
bility still give significant deviations at shock 
pressure ratios near 10. The present results and 
the results of reference 2 are in agreement within 
less than 15 percent for shock pressure ratios up to 
10 and then diverge to a 20-percent variation at 
a shock pressure ratio of 20. 

Figures 5 and 6 show that the cold gas con- 
tributes only a small part of the total attenuation 
for /WPoo<20. At Pdo!p^ = 20 the relative a 
contribution is larger than that for p& u /p <x> < C2() ] 
however, it is only about 4 percent, 15 percent, 
and 25 percent of the total for m— 1, 7, and 13, 
respectively. 


EVALUATION OF NONLINEAR ATTENUATION EXPRESSIONS 


The expressions derived previously are based 
on a linearized or small perturbation analysis. 
However, for many conditions encountered in the 
shock tube, the attenuation is no longer small. 
In order to maximize the total available experi- 
mental testing time, most experimental work is 
done at values of l nearly equal to the total length 
of the low-pressure side of the shock tube. At 
these large values of / the shock strength often 
has decayed markedly from its value for small /. 
Consequently, relations for the attenuation under 
these conditions would be very desirable. 

An approximate method to obtain the attenua- 
tion for the cases where the small perturbation 
analysis is invalid will be described. First, con- 
sider parameters P and Q which are related to P 
and Q by 


ylt 


(43) 


(The parameters P and Qof tins report are identical 
to the parameters P r and Q' of ref. 1.) When 
equation (43) is substituted info equation (60) 


of reference 1, the following equations result: 


/*■ A -S 

a t dt 

e 5 1 J 


_1 a 1 bp 1 bU 
~ya i p bt dt 


a, 8t 

1 dQ 


=|^ 7 o / { i( 7-1)A/T1 ] A/ 

T — T 

j — 2/3 1 w 1 aw 


+ N Pr ~ w - T 


} 


i di‘ 


<Pj> 

a, dt 

II 

fel- 

1 8Q 
J>‘ 5t j 

06 

V J 


(44a) 


(44b) 


(44c) 


As discussed in reference 1, the changes in P and 
Q are evident in changes by wave motion of veloc- 
ity and pressure but do not indicate the various 
changes in entropy. For example, if equation 
(44a) is solved for bpjp } 


dp_ya, ( sP | ( 45 ) 

p 2 a \« e «. / 

A 

The value of P is associated with waves moving 

A 

with the flow at a velocity of u-\ a whereas Q is 
associated with waves of t lie opposite family mov- 
ing at a velocity u—a. Now, for the linear at- 
tenuation theory it is assumed that reflections at 
the shock wave and the deviation in entropy rise 
across the shock wave may be ignored; that is, 

A A 

Q 9 * — Qv»—Q&o or bQ,;s=bQ vs -= 0. Consequently, an 
alternate form for the attenuation expression is 


&P_ P*>—Peo = y ( h I" ( ?>P\ , ( &P\ , ( 

P Ppo 2 (l [_V ) a \ & t / reflection \ ^ 

(40) 

where the three terms on the right-hand side oT 
the equation represent, respectively, the contri- 

A 

but ions of region a, of the Q wave generated in 
region ft and reflected at the entropy disconti- 
nuity, and of the P wave generated in region 
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Consider now the last term only. The incremen- 
tal change bP along the characteristic is then 



a. bt dl ~T) 



L Cf4>p t $bt 


(47) 


A change of variable may be made (eq. (21) of 
ref. 1) to £ since 


Ubt 




( bP\ r««.> 

\ Cl e / p Pj € (tf) cr J\ f s A'/jj 1 


Cf4>r,pd£ 


f 'bP\ \ Cta vM s —M& 

\Cl e )p T)J f) <J M s Afp 1 


4>r,P c f,p(£)d£ 


(48a.) 



aM s —Mp / 
aMg—Mp— 1 \ I> / 


IV/ ,»<!% 

(48b) 


Substituting equation (45) into equation (48) 
yields: 


' Pr»—P ea\ /«A 
, PPo )p,r V P )p,p 


(49) 


In equation (49), the left-hand side represents 
the pressure perturbation at the shock wave due 
to the wave generation only along the forward 
running (slope of u-\-a) characteristic. For a com- 
plete linearized treatment K and may be taken 
outside the integral and equation (35b) may be 
employed t o obtain : 



Thus, K is the ratio of the contribution of the 
waves generated along the forward running charac- 
teristics to the total waves generated in region /?. 
The remaining portion of the attenuation contribu- 

A 

lion of j 8 results from Q waves reflected at the 
entropy discontinuity and will be designated by 


| ) which is equal to 

\ P /»A+s 

(?) . =(1 -X)(^=^f) (51) 

\P/0,Q+3 \ PPo /P 

The value of K is plotted against shock pressure 
ratio in figure 7. The fact that K does not depart 
significantly from 1.0 means physically that the 
principal contribution to attenuation in region (3 

A 

arises from the P waves. Since figures 5 and 6 
have shown region 0 to have a much larger effect 
on attenuation in general than region a, it is obvi- 
ous that the theoretical dominating factors for 

A 

attenuation are the P waves of region /3. This con- 
clusion has been discussed previously in references 
1 and 2. 

A 

Since the F waves of region (3 dominate the lin- 
earized attenuation solution, it is next assumed 



Future 7. Ratio of linear attenuation effects generated 
along forward running characteristic to total attenua- 
tion effects of region as a function of shock pressure 
ratio. 
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that a correction for the linearized solution may 

A 

be found by operating only on the P waves in 
region 0. Thus P waves generated in region a and 
transmitted at the entropy discontinuity as well 
as the Q waves generated in region 0 and reflected 
at the entropy discontinuity will retain their origi- 
nal linear or small perturbation values even though 
the attenuation is no longer small. It is further 
assumed that region 0 may be subdivided into a 
number of smaller regions in each of which the 
linear attenuation relations for P are valid. This 
treatment is illustrated in figure 8. The arbitrary 
interval A l determines the x-wise extent along the 
shock wave of each of the regions designated (T), 
.... Each of these regions is bounded by the 
shock wave and two fluid particle paths where each 
fluid particle velocity is equal, respectively, to 
that generated by the shock wave at the beginning 
of each new interval. The inviscid flow inside 
each of these regions is considered to be constant ; 
and, consequently, there is a small discontinuity 
in the inviscid flow across the particle-path bound- 
aries assumed in the model. These discontinuities 
can not, of course, exist in the actual physical flow 
which requires a continuous variation throughout 
all the regions as well as reflections from the shock 
wave. The errors introduced by the assumption 
of constant quantities in each region are not con- 
sidered to he large and should be of approximately 
the same order as those found in the familiar 



Ficfre 8 . — General wave diagram for nonlinear analysis 
in region p. The short-dashed line is an extrapolation 
of the shock-wave path and the dotted lines arc extrap- 
olations of the characteristics. 


steady-flow graphical characteristic solutions of 
finite mesh size. 

In order to simplify the computational proce- 
dure, it is assumed that for a given region the 

A 

slope 8x/8t of the P characteristic and of the shock 
wave are constant hoth inside and outside of that 

A 

region. Tims, in figure 8, when the P contribution 
of region ® to attenuation between A / and 2 A / is 
computed, the assumption is made tha t the charac- 
teristic 6,7 and the shock path 0,1 may be ex- 
tended to intersect at point 2. (Numbers refer to 
points in fig. 8.) The corresponding corre ct re- 
gional characteristic line and shock paths are 6,7,8 
and 0,1,8 which are shown in this illustration as 
intersecting also at the same value of x as 2. This 
intersection at the same value of x is only an ideal- 
ization and is not the true physical picture in gen- 
eral. However, since the attenuation effect (gen- 

A 

oration of P waves) falls off rapidly with distance 
behind the shock (similar to the fall off in local 
skin friction with distance back of a sharp leading 
edge in steady flow), the contribution to attenua- 
tion in the interval from A / to 2 A / due to generation 
along 6, 7 is much less than that due to genera tion 
along 7, 8; thus, small errors in the location, of 6,7 
will result in very small errors in the attenuation 
at 2A/. This assumption for establishing the inter- 
section points of the characteristics and the shock 
wave downstream from a region (§) without knowl- 
edge of the downstream shock-wave attenuation 
permits the easy computation of the influence of 
region (n) for all downstream shock locations. 

When the regional approach described above 
is applied, the attenuation for the first interval 
Al is identical to the complete linear approach. 
Thereafter, however, the various second-order 
effects are felt. The effect of region ® on attenu- 
ation of the shock during the interval from Al to 
2 Al differs from its effect in the basic linear 
theory because of the convergence of the particle 
paths since U u <CJh* This may be shown as 
follows. From equation (48), 


6P\ | 7 = J_ Hr/ aM-My 
. Pjh \oA/* A 'I p 



dZ 


(52) 


The numbers refer to the points on fi gure 9 (a). 
In this figure the lines 5, 1, 10, 11, and 6, 7, 2 are 
drawn with slope (fj+ah; the li nes 7,10 and 
0,9, 1 , 11,2 with slope T7 sl ; the line s 0, 5, 10, 6 and 
977 with slope Ui] and the line 1,7 with slope U u . 
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(a) Details for computations of region I. 

(b) Details for transitional flow. 

Fiotre 9. — Detailed wave diagrams for nonlinear analysis 
in region p. Short-dashed lines indicate extrapolations. 


Therefore, £ 7 =9, 7=0, 10 = £,„. Thus, 

^SP\ 7 _ 1 ft../ aM-Me \ 

V «. A e Dj lA > ^ 

_ 1 / \ / ft« 

nW-M^-iAUo 

(53) 

Thus, the P contribution of region ® to attenu- 
ation between 1 and 2 is (from eqs. (45) and (53)) 

d)„. 

-(&i) m 

\p^/l.0tol a \pf3 / 1, o to in V 

Tlie substitution of equation (40), with g„ replaced 
by g n , into equation (54) yields 

2 _ n +1 

^ \ /«oo^\ "+3// i \„+3 


Af-G-f] 

^..[Gf-Gf] 


(55) 


where f — ) is the linear attenuation'bf a'shock 

VPco/lAl 


A 

due to P effects over the initial interval Al in region ©. 

The relationship between l n and l x is derived in appendix B and is 


Consequently, 


7 T Z '(H 

(fl „ „ G-'f*] ] 

/sp\ j sj\ f/^_ 

\i>Ji, h to , 3 { i) \h J J 


(56) 


(57a) 
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which for equal intervals A/ becomes 



I. o to 


I. hi to 


=( -) 1 

- n+i »±n " 

J2i‘+3 — (Zj)«+sJ 

/ £ 7)\ 

r 71+1 «+jn 


[_drt+3 — (2Zi) »+ ; *J 

f 8 \Pa>/l,M 

> 


( 571 >) 


In this form the nonlinear attenuations can easily be computed. 

In this manner the influence of region (f) on attenuation at any desired value of l ma^ compute t 
once (lie attenuation at U has been found. The influence of region © is found in a similar manner 

by shifting the effective origin of the coordinate system to h and finding (^) n for tlu ‘ attenuated shock 


strength at l\\ that is, 


n+ 1 
r+3 


Xpo/llJx to l 2 \poo/U A l \ V a> / \ iJ / 

m -m 

\p<x> 'll.U to *3 V/^/lI.A! Wk—ll/ \h—ll / j I 


J 


or for equal intervals of A l, 


( S P) =(—) [2 J+»—(Z, ,)"+*] ] 

\p<x>/n.i 2 to ^ \P<x>/iiAi I 

(OPS) =(*£. \ [:t«S-(2Z 1I )^ i ] 

\/>=./iu 3 to i 4 \p co/ii,A( 


(58a) 


(58b) 


(5 Sc) 


The total nonlinear attenuation at a distance l from the diaphragm station which lias been sub- 
divided into several intervals &l is then expressed as the sum of the linear contributions of region a and 
the reflections from region ft added to the nonlinear contributions of region ft. T he following expression 

is obtained for the nonlinear attenuation: 



(59a) 


(59b) 


where i is the index. 
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Second-order attenuation in the presence of tran- 
sition may also be treated by this regional system. 
The contributions of a and of Q and S in are still 
treated as entirely linear but including transition. 
Equation (39a) applies in its entirely to the a con- 
tribution; and, when the first term on the right- 
hand side is multiplied by \—K m and the second 
term by 1 — K n , the resulting equation gives the 
Q and S contribution of region /?. 

The transitional 7 J contribution is treated as 
follows: Let l* be the position of the shock at 
which the transition in the flow behind the shock 
first affects the shock. The wave diagram for 
transition is shown in figure 9(b) where for il- 
lustrative purposes it is assumed that /*— 2 A/. 
Xow ? in order to retain the facility of computation 
afford ed by the regional system with equal At and 
a constant l*, the value of It* will, as a result, 
vary slightly from region to region. Since 




the magnitude of this variation 


may be found by examination of figure 10 which 


shows the parameter 


W Z*/v 


plotted against shock pressure ratio. From this 
figure it is evident that, if the shock pressure 
ratio should attenuate, for example, from 20 

to 15, or from 5 to 4^ there would be about a 

10-percent decrease in It* for a given l* and 
The errors introduced by such a varia- 
tion in It* are not deemed to be important enough 
to force the abandonment of the equal At com- 
puting scheme. For the remainder of this paper, 
It* will be taken as the value of the transitional 
Reynolds number in region (T ). The bp contribu- 
tions for the transitional case are expressed by 
the following equations which are modifications 
of equations (57), (58), and (59) (the subscripts 


m and n 

refer to 

conditions before and 

after 

transition, 

respective 

ply): 






2 m-f 1 


m 

\Pa oAv, m,M 

({7/«^/«)jv 

/ 

V V 00 / 

v m+ 3 /A/\m +3 

' \d) 

(60a) 




■ 2 n-f -1 


m. , 

\P oo / A , n, A l 

= (f)AXv i 

Cf) 

n+ 3 /A /\«+3 

1 \d) 

(60b) 



Fku re 10. — Regional transition Reynolds number per 
foot of shock-wave travel per atmosphere of ambient 
pressure plotted against shock pressure ratio. Circled 
symbols denote regions. 
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h> N ~'+s 

NUMERICAL EVALUATION OF NONLINEAR THEORY 

Several computations to determine the non- 
linear correction factors for attenuation were per- 
formed for values of the interval Al/D of 3.5 and 
14. These particular values were chosen because 
they represent increments of Al of 0.5 foot and 
2 feet for the shock tube employed in the experi- 
ments to be described in a later section. Typical 
curves resulting from such computations are 
shown in figure 11. The ratio of the nonlinear 
attenuation to the linear attenuation is plotted 
against l/D for an initial shock pressure ratio of 4 

and a value of of 0.01 X10 6 . For values of 

^co 

l/D^>50 the curves for ~=3.5 and 14 agree within 

about 1 percent for the laminar case and within 
about 3 percent for the turbulent case. The low 
value of a^D/v^ accentuates any variations 
between the two computations; thus the case 
illustrated gives a discrepancy near a maximum 
rather than near a minimum. Examination of 
several such pairs of curves resulted in the con- 
clusion that the slight increase in accuracy 

obtained by using ^=3.5 did not justify the 



(a) Turbulent flow in «; laminar flow in 0. 

Figure 11. — Ratio of nonlinear to linear attenuation 

plotted against shock-wave position for ~°=4.0 and 

P oo 

a D 

-^-=0.01 X10 B . 



(b) Turbulent flow in both a and fi. 
Figure 11. — Concluded 


fourfold increase in labor. Consequently, the 
computations with ^=14 are used to predict the 

nonlinear at tenuation for l/D*? 50. 

It is obvious that the finite size of Al will 
introduce errors in the ratios NL m%n jL mt n which 
are largest near l- >0 since the nonlinear and 
linear attenuations are identical for the first 
interval. (See fig. 11.) However, the errors 
introduced in the attenuations NL m n themselves 
are small since L m>n -> 0 as l— >0. To represent the 
physical flow in this region accurately would 
require that Al approach 0. Interpolation for- 
mulas giving acceptable accuracy near are 

assumed to have the form 


Nonlinear attenuation fr ^ , ,\ 

— — : = (Constant) 

Linear attenuation 

Nonlinear attenuation _ / ^ o ^^^ 
Linear attenuation 



Y (62) 


for ^0<2^5G^ for laminar and turbulent flows, 

respectively, since the linear attenuation is 
proportional to (l/D) 0 5 and ( l/D) 0 S for the laminar 
and turbulent flows. 

The constants are chosen to match the computed 
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curves for ^—14 at -^=56. From figure 1 1 it is 

evident that the errors resulting from the appliea" 
lion of tins interpolation formula are less than 

the aforementioned errors at ^=5 0 and are 

thus acceptable. 

In order to obtain curves of the transitional 
nonlinear to linear attenuation ratio NL 7iT /L 7tT 
for a constant value of /?* and several shock 
pressure ratios, a cross-plotting technique was 
used. At a given shock pressure ratio the values 
of NL 7tT lLj T were computed for several values of 
the ratio l*JD to Al/D. Since each /* represented 
an /?*, the resulting ratios for each IjD could be 
plotted against R* for a given shock pressure 
ratio. The values for a particular /?* could then 
be read from these plots to produce a master plot 
with a common value of /?*. 

PI o t s o f t h e r a t io o f nonl i n ear a 1 1 enu a t i o n t o 1 i n ear 
attenuation are shown in figures 12 to 15. In 



(a) ^=25. 


(c) £=75. 

Figure 12. — Ratio of nonlinear attenuation to linear 
attenuation for region a turbulent and region laminar 
plotted against shock pressure ratio. The shock wave 
is located at a distance IjD from the diaphragm station. 

5G97GS — Cl 4 



Figure 12.- Concluded. 


these curves region a is always considered as having 
turbulent flow whereas four cases are considered 
for region /?: namely, (a) laminar flow, (b) turbu- 
lent flow, (c) transition with /?*= 1.25X 1 0 6 , and 
(d) transition with 7?*=2.5X10 6 . The cross- 
plotting parameters are shock-tube Reynolds 
number a^Djv^ and length of shock-wave travel 
expressed in hydraulic diameters IjD. At the 
lower pressure ratios p&o!p«>, curves for more 
values of <i m Djv„ are shown than at the higher 
pressure ratios. This limitation resulted from the 
considerations of the restriction of the validity 
of the theory to an ideal gas, the region of experi- 
mental data of this report, the most likely general 
region of experiments for other facilities, and 
priority for computing effort. Since figures 12 to 
15 arc the result of cross plotting, the accuracy is 
assumed to be about 2 percent. 
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(b) 2 ;= 50. 


Figure 13. Ratio of nonlinear attenuation to linear 
attenuation for both regions a. and p turbulent plotted 
against shock pressure ratio. The shock wave is lo- 
cated at a distance IjD from the diaphragm station. 

An analytic closed-form investigation lias been 
made of the fact that the limit of NL 7t ijL 7t \ 
approaches 0 when p$ a !p<x> approaches 1 whereas 
the limit of NL 7J /L 7J approaches 1.0 when 
PuJPn approaches 1.0. This second -order analytic 
solution for weak shocks indicates that a value 
of n=3 in /3 is a critical value; all solutions with 
n>3 approach a limit of 1 and those with n<3 
approach a limit of 0. Of course, in all eases 
the absolute value of both the linear and nonlinear 
attenuation must approach zero as ppolpa, ap- 
proaches 1. Since a has only a secondary effect 
on attenuation and since the expansion fan has 
been replaced by a “negative shock/' the refine- 
ment of transition in region a was not deemed 
necessary. 

RESULTS AND DISCUSSION 

EXPERIMENTAL APPARATUS AND PROCEDURE 

Extensive shock-attenuation data were obtained 
in a high-pressure shock tube 2 inches high by 



tc) 1=125. 

(f) 1= 150. 

Figure 13. Concluded. 

1 1 '2 inches wide in the Langley gas dynamics 
laboratory. This shock tube is the same as that 
described in reference 1 with added velocity- 
measuring equipment. Light screens placed at 
eight stations in the low-pressure section made 
possible measurements of the complete distance- 
time history of the motion of the shock wave for 
a wide range of shock pressure ratios and flow 
Reynolds numbers. Figure 1G shows schematically 
the arrangement of the shock tube and associated 
equipment, and the low-pressure section with 
associated optical systems is illustrated in figure 17. 
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Figure 14. — Ratio of nonlinear attenuation to linear 
attenuation for region « turbulent and region j3 
transitional plotted against shock pressure ratio. 7?* 
= 1.25X10 fl . The shock wave is located at a distance 
JjD from the diaphragm station. 

Air at room temperature was used in both high- 
and low-pressure sections. The normal arrange- 
ment for high-pressure section air supply and 
low-pressure section vacuum systems is shown in 
figure 16. In a limited number of low-density 
runs, evacuating the high-pressure section was 
necessary; for these cases, an auxiliary vacuum 
system, identical with the normal low-pressure 
system, was substituted for the pressure system 
shown. Conversely, certain high-density runs 
required pressurization of the low-pressure section; 
in this case a simple pressure system replaced the 
normal vacuum system. The pressures in both 
high- and low-pressure sections were adjusted for 
each run. All data were obtained from tests 
where the diaphragm was punctured by a hand- 
operated plunger. In t his way diaphragm pressure 
ratio and consequently, theoretical shock-pressure 
ratio Ppo/Poo were closely controlled. Bourdon- 
lube gages were employed for pressure measure- 
ments, and the vacuum systems utilized a modified 



(f) j-=l50. 


Figure 14.- Concluded. 

barometer for pressures in the range from 4 inches 
mercury absolute to atmospheric pressure, and a 
0 to 100 millimeters mercury absolute pressure 
gage was used for the low pressures. 

Static pressures were structurally limited to 
1,000 pounds per square inch gage in the high- 
pressure sections and 250 pounds per square inch 
gage in the low-pressure sections. A vacuum 
limitation of about 0.01 atmosphere absolute 
pressure resulted for the low-pressure section be- 
cause the light screen systems became unresponsive 
for (lie low pressures. In the high-pressure section, 
low pressures were limited by failure of dia- 
phragms to burst properly. 

Diaphragms made of thin metal foil were used 
for the low-pressure runs. The most useful ma- 
terials were soft aluminum foil, 0.001 inch thick, 
and soft brass foil with a nominal thickness of 
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Fictre 15, Ratio of nonlinear attenuation to linear 
attenuation for region « turbulent and region /? 
transitional plotted against shock pressure ratio. R* 
= 2.50X10*. The shock wave is located at a distance 
IjD from the diaphragm station. 

either 0.00125 or 0.0015 inch. Willi these ma- 
terials, it was possible to obtain good bursts for 
pressure differences across the diaphragm ranging 
from 10 pounds per square inch to 100 pounds per 
square inch. Good bursts for the range of pressure 
difference from 100 to 1,000 pounds per square 
inch wore obtained by using spring-tempered brass 
shim stock with thicknesses ranging from 0.008 to 
0.021 inch and scribed to various depths in an 
x -shaped pattern along the diagonals of a rectangle 
representing the shock-tube cross section. When 
punctured under pressure, the diaphragm split 
along the scribe marks, and the four triangular 
pieces of material folded back against the wall and 
presented minimum resistance to the flow. All 
other conditions being equal, runs where this type 
of burst took place resulted in minimum shock 
attenuation compared with attenuations resulting 
from bursts where metal or acetate-type dia- 
phragms wore shattered. Material thicknesses 



(f) 7)= 150. 


Figure 15. — Concluded. 

and scribe depths were determined so that punc- 
turing pressure was just below the pressure at 
which the diaphragm would have failed. The 
unscribed foil diaphragms used for small pressure 
differences split along the diagonals in this same 
way when punctured at the center, and optimum 
results were also generally obtained with these 
diaphragms. 

Shock velocity measurements were made with 
miniature sehliorcn systems located at eight posi- 
tions in the low-pressure section. Figure 16 
shows schematically the position of these systems. 
The optical and electronic systems were essentially 
the same as those used in reference 1, in which 
the signal generated in a photomultiplier tube by 
deflection of a beam of light upon the tube was 






-93040 

Figure 17. — Photograph of experimental setup showing shock-tube low-pressure section with velocity-measuring stations. 


amplified and used to trigger a thyratron. The 
thyratron output pulse, in turn, started or stopped 
a counter Chronograph. Figure IS illustrates 
one complete optical system, including the chassis 
containing the photomultiplier-amplifier-thyratron 


circuit, which is shown in figure 19. Wooden 
shields were employed to keep stray room light 
from falling upon the photomultiplier tube. 

For the multiple systems employed herein, each 
thyratron output pulse was channeled to two 
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Figure 18. — Photograph showing details of velocity-measuring optical system. 


chronographs. These chronographs indicated 
si lock traversal time between any two adjacent 
stations, and any one of three types of chrono- 
graphs, 8 megacycles, 1.6 megacycles, and 1 mega- 
cycle, was employed for each position. The over- 
all time interval between the first and last stations 
was measured with a 100-kilocycle chronograph 
for a check upon the sum of the individual measure- 
ments. 

REDUCTION OF EXPERIMENTAL DATA 

The experimental shock-wave pressure ratio 
was computed from the measured time interval 
At for the shock to pass between two measuring 
stations a distance Al apart from the relation 


p** = 27 SAl/At\ 2 7—1 
pc 7+1 \ / 7+1 


(63) 


Tli is value was assumed to represent the shock 
strength at a position midway between the two 
stations, and the maximum error was estimated to 
be less than 1 percent. The theoretical shock 
pressure ratio p& 0 !p<x> was computed from the 


diaphragm pressure ratio just prior to burst. For 
an ideal gas the maximum error in pp 0 was esti- 
mated at 0.1p OT (for a shock pressure ratio of 10 at 
p„ = 0.005j? ro(5W ); and the maximum deviation in 
the ratio of computed (p&/Pc)o to true (p&/p m )o 
was estimated to be approximately 1 percent. 

COMPARISON OF THEORY AND EXPERIMENT 

The general method of comparison between 
theory and experiment will be to compare the 
measured and predicted attenuation for particular 
values of p&olPc and a m Djy „ on individual curves. 
However, it is of interest to first consider a few 
typical curves where the data for a constant value 
of pno/pa, but with varying values of a^Djv^ are 
shown on a single plot. Such plots are shown in 
figure 20 for values of p$ 0 }pc of 4.0 and 10.0. The 
experimental data are averages of several runs on 
a given day and the number of runs for each data 
point is indicated on the figure. On some runs in 
which the density change across the shock wave 
was small, all the velocity-measuring stations did 
not register because of variations in their sensi- 
tivity; and, as a result, there are gaps in the 
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Figure 19. Schematic diagram of photomultiplier and output stages for shock velocity-measuring system. 


experimental data. ^For example, see fig. 20(b) 

at (t ^=0.005.^ 
v*, / 

One fact immediately evident is the nonre- 
peatability of some of t lie data, even when com- 
pared on a daily average basis. An allowance 
for an error in pa/p^ of ± 1 percent in the experi- 
mental and data-reduetion technique will bracket 
most of the observed discrepancy in the averages; 
but certain runs at very low values of <i w T)jv „ still 
fall outside this range. 

If an attempt were made to extrapolate a curve 
from the data points to Z=0, an inflection would 
often be required in the curve between 1 = 0 and 
Z=G to make it pass through the theoretical value 
of P&o/Pn at l— 0. In order to illustrate this point, 
connecting lines have been drawn in figure 20 for 
some of the values of Similar behavior 

is found in the experimental data reported in 
figures 11 to 14 of reference 8. Now all the at- 


tenuation theories based on wall effects which arc 

ip 

known to the authors predict UJ ^>0 in regions 

of laminar or turbulent flow. At the transition 
point the theory of the present paper usually 

pmli.-ur-'fcfc) >(%'"-) <0. 

\ m / laminar \ / transition 

Consequently, if this inflection is to arise from 
wall effects, if must be caused by transition. The 
interferometric studies of reference 1G have found 
values of /?* of about 1.4X10* and 2.0X1 0 6 for 
zWp«> — 2 and 24, respectively. Since values of 
77* of the order of 0.5X10 6 are required to cause 
an inflection at Z=G for p,j < ,/p a> = 4 and 10, it docs 
not appear likely that transition is the cause of this 
inflection. This inflection will be discussed more 
fully later in this paper. 

The unexpected inversion of Reynolds number 
effects for the lower values of l should be noted. 
For example, at /= 5. 2 and p^ 0 jp m — 10 the attenu- 


at ion for 


=0.005X10° is about one-fourth that 



Experimental local shock pressure ratio 
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Figure 20.— Concluded. 


the predicted attenuation when the attenuation is large but there is also a marked decrease in the 
predicted Reynolds number effect on attenuation. The following table based on figures 21 and 22 for 
1—22 feet illustrates these facts: 


Curve 

Figure 

van.es or j£-« 

for values of ^ of - 

v co 

Values of (—Vffr 1 ), for ^-10 

\P <x>/ VPoo// = 22 Poo 

for values of a of — 

^co 

R * 

0.01 XIO 8 

o.ixio* 

1X10 8 

o.otxio* 

0.1x10* 

lXIOfl 

Lt.x 

21 

1.45 

0. IS 

0.17 

2.62 

0.66 

0. 10 

00 

U. 7 

21 

2.30 

1.45 

.91 

6.10 

3. 85 

2.43 

0 

NX7.I 

21 

1.09 

.44 

. 16 

1.87 

.58 

.10 


NTa; 7 

21 

1.39 

1.04 

.74 

3.22 

2.50 

1.83 

0 

AT.-.r 

22 

1.09 

,99 

.73 

2.65 

2.46 

1.8 

1.25X10* 

NL-.t 

22 

1.09 

.92 

.72 

1.87 

2.38 

1.8 

2,5X10* 


The nonlinear attenuation predictions for values of Ii*= 1.25X10® and 2.5X10* are shown in 
figures 22(a) and 22(b). Since, for these values of /?* Z*=0.122 and 0.244 foot for ^^=10* V&o 

F«> Pco 

= 10, the nonlinear curve for turbulent flow is nearly identical to that for transition and the turbulent 
curve is used in tin* figure. These curves and the table indicate that transition at a constant 
value of R* appears to decrease further the spread of attenuation with a^D/v^. The discon- 



32 


TECHNICAL REPORT R-85— NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 



Fi(;rRE 21. ---Theoretically predicted local shock pressure ratios plotted against shock distance from diaphragm. 


tinuous change in d ^ 1 ^jj 00 at Z* is obvious. An- 
other point of interest is the predicted variation 
of NL Jt r at 1= 22, pvofp a = 10, and 7?* = 2.5X10 6 ; 
the attenuation is less for of 0.01 X10 6 and 

1X10 6 than it is for a^D/v^ of 0.1X10 R . The 

completely laminar flow of — — O.OIXIO 6 and 

the nearly completely turbulent flow ol ~ — = 
lX10 fi result in almost identical attenuations. 
Tlio flow of —^—=0.1 X 10 fi is mainly turbulent 

V ^ 

(Z* — 2.44) and, as expected, the attenuation is 

greater than for ^-=1 X10 6 . Thus, in this case, 

transition has resulted in an inversion of attenua- 
tion with Reynolds number under a certain set of 

conditions. However, for 7?*=1.25X10 6 ,^=10, 

and Z— 22 feet, the pattern of increasing attenua- 
tion with decreasing a^Djv^ is once again evident. 

In figures 23 to 36 the variation of the average 
experimental shock pressure ratio with distance 
is compared with the various attenuation predic- 
tions of the theory presented earlier in this paper. 


The nominal values of Pf*„lp w range from 1.5 to 
10 and aj)!v m covers the range of 0.005X10 6 to 
15X10 6 . The hydraulic diameter of the shock 
tube is }'7 foot and the maximum value of l at 
which Pcs/pc o could be determined was 17%> feet; 

thus a maximum experimental value of 


result cal. 

The marked improvement obtained from use of 
the nonlinear theory when the predicted attenu- 
ation is large becomes evident upon inspection of 
figures 23 to 36. The unusual behavior near Z->0 
which was mentioned earlier can now be examined 
more closely. For the higher values of shock 
pressure ratio (8, 9, and 10) and low values of 

^^^0.01 X10 R , the measured attenuations at the 

v ^ 

first station are much smaller than any of the 
theoretical predictions of this report or of refer- 
ence 2. In fact, for ^-^-=9.0 and — — 0.005 XI 0 R , 

p oo r o 0 

one set of averages gives a negative attenuation 


( \ of such size that even the estimated 
\P» Pa>/ 

1 -percent -error margin is not sufficient to make 
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Figure 21. — Concluded. 


at the first station (^=36.()\ In 

P~ P~ \D ) 

reference 9 it was also found that for high values 
of P&olPn the maximum shock velocities were 
greater than the velocity theoretically computed 
for an in viscid fluid both with and without con- 
sideration of variable specific heat and gaseous 
imperfections. Hydrogen and helium were used 
as the driver gas and air as the low-pressure fluid. 
The maximum shock velocity occurred at about 
40 to 50 diameters from the diaphragm. 

On the other hand, references 4 and 8 did not 


find any values of — for values of p$ 0 !p<x> in 
P ca P oo 

the same range as the present experiments. 
These works proposed a “formation decrement” 
defined as “the difference between the Rankinc- 
Ilugoniot shock strength and the maximum shock 
strength obtained after the formation distance.” 
(See ref. 8, page 17.) This decrement was then 
attributed to the imperfect diaphragm burst pro- 
ducing a series of compression waves which 
eventually coalesce to form a si lock weaker than 
that for the case of a theoretically perfect burst. 
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Distance from diaphragm, ft 


Figure 30. — Comparison of averaged experimental local shock pressure ratios with various 
theoretical predictions. ^=4.0. The vertical lines in (a) and (b) indicate the range of 

scatter in an averaged group when the scatter exceeds the symbol size. Flagged and un flagged 
symbols are used to distinguish between averages from different days. 
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For the lowest shock pressure ratios at the 
highest shock tube Reynolds numbers, the opposite 
trend (measured attenuation larger than theory) 
is apparent. An extreme example of this is the 
attenuation measured at the first station for 

— = 15X10® and ^=2. (See fig. 25.) Three 
V <x> P 00 

of the four runs used in obtaining this average 
value showed this behavior, which is attributed 
to the formation process. 

The processes giving rise to the behavior for 
1/D-+0 are not understood. The experimental 
data of this report generally show a trend of de- 
creasing pt a fpfi 0 with increasing a^Djv^ at the first 
station. If this effect were to be explained on a 
viscous basis, it would appear that an opposite 
trend would appear; hence, the diaphragm burst 
is probably the governing factor. The bursting 
phenomena are in turn governed by the diaphragm 
materials (weights and rupture characteristics) as 
well as the pressure load. The diaphragm opening 
time (time for the diaphragm sectors to fold 
against the wall) was estimated by a method 
which has previously been checked experimen- 
tally. No correlation was found between the 
opening time and the behavior near /— >0. It was 
noticed, however, that, in the experimental runs 
which exhibited the most marked inflections, 
foil diaphragms were used. On the other hand, 
some runs exhibiting very little or no inflection 
also used foil diaphragms. Therefore, it does 
not appear that foil diaphragms are solely re- 
sponsible for tlie inflected data points. 

Regardless of whether the maximum experi- 
mental value of prs/pn is greater than or less 
than the ideal value, the behavior near the dia- 
phragm station is not governed by wall effects 
but by the diaphragm burst and the resultant 
three-dimensional flow first established. This 
highly rotational viscous flow does eventually 
become essentially two-dimensional, with the 
exception of mixing and vorticity at the interface 
between the driver and driven gases. Conse- 
quently, any attenuation theory based on wall 
effects cannot predict the initial behavior near 
the diaphragm station. As the distance from 
the diaphragm station increases, the ratio of tlio 
influence of the initial bursting flow to the influ- 
ence of the integrated wall effects decreases; 
therefore, the physical variation of p c Jp m should 
approach the theoretically predicted attenuated 
value asymptotically as l/D increases. 


The experimental data of figures 23 to 3G 
approach the nonlinear theory (considering transi- 
tion) in a manner very similar to that just de- 
scribed. For high values of p&olp m and low 
values of a m D/v m where the effect of the forma- 
tion process results in a significant inflection, the 
deviation from the nonlinear curve persists to 


the larger values of l/D, For 10, 9, and 8 

with —"—=0.005 X 10 6 , the formation effect was 
so large that the experimental data never ap- 
proached closely the nonlinear curve for -4^ 125. 


As a^D/vn increased and the formation effect 
decreased, agreement between experiment and 
nonlinear theory improved both in convergence 
of the experimental and theoretical values at a 
lower value of l/D and in maximum deviation at 
the highest values of l/D . 

The nonlinear curve for 7?*=2.5X10 fi appears 
to agree more favorably with the data over most 
of the range of p$ 0 !pcx>* At the lower values of 
p fio! P co there appears to be a tendency for the data 
to depart from the J?* = 2.5X10 6 curve and ap- 
pro ach the 7?*=1.25X10 6 curve. The inter- 
ferometric measurements of reference 16 indicate 
such a trend of increasing R* with increasing 
P$o!p<»* This trend might also be expected from 
comparison with steady-flow experiments since the 
wall cooling increases as p&o/p* increases. 

The comparisons of figures 23 to 36 between the 
nonlinear transitional theory and the experimental 
data show that this theory is valid for the pre- 
diction of experimental attenuation except for the 
lower shock tube Reynolds numbers at the higher 
pressure ratios. These latter conditions are those 
under which it appears that the shock-formation 
processes dominate the entire flow. 

As stated previously, errors of unknown mag- 
nitude were introduced by the averaging of wall 
effects across the flow. The fact that the present 
theory was able to predict fairly well the measured 
attenuations over the entire ran^e of shock pressure 
ratios and Reynolds numbers (except near the 
diaphragm) indicates that the errors introduced 
in the averaging process are either not serious or 
else self-compensating. 

The range of boundary -layer thicknesses in 
region 0 (which has the predominant influence on 
attenuation) was determined by methods similar 
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to those of references 14 and 15 for laminar and 
turbulent flows, respectively. Tn the experiments 
reported heroin, with the shock 20 feet from the 
diaphragm, the maximum bound ary -layer thick- 
ness varied approximately from O.OOG to 0.25 
inch for laminar flow and from 0.07 io 1 inch for 
turbulent flow. When compared with the 0.75- 
inch half-width of the shock tube, the turbulent 
thickness was significant for much of the experi- 
mental data. On the other hand, laminar bound- 
ary-layer thickness was generally small. Tlius 
the averaging process as used for this theory 
appears to apply equally well to thick as well as to 
thin boundary layers. 

MAGNITUDE OF OTHER NEGLECTED SECOND-ORDER 
EFFECTS 

Since the treatment presented lias considered 
only one source of nonlinearity, it is advisable to 
examine briefly the influence of other neglected 
sources. The three most obvious factors neglected 
are the finite extent of the expansion fan, the 
nonlinear effect of the wave interactions at the 
entropy discontinuity, and the reflections at the 
shock wave. 

The effect upon attenuation of treating the 
expansion wave as one of finite width with varying 
free-stream properties was calculated by using 
the results of reference 14 for pressure ratios 
pJ Pco of 1.0, 2.9, and 4.5 (£,« of -0.G, 0, and 0.5, 
respectively, in the notation of ref. 14). Com- 
putations were not made for higher shock pressure 
ratios because the solutions in reference 14 did 
not extend above £ fe =0.G. The laminar-flow 
wall-shear and heat -transfer distributions through 
the cold-gas regions given by this reference were 
used to compute the skin-f fiction coefficient 
which was, in turn, utilized to compute the 
attenuation. 

When the method of the present paper was 
used, it was found that, for the three cases com- 
puted, the net change in attenuation through 
approximation of the finite expansion by the 
negative shock was practically zero for ppo!p<x> 
equal to 1.6, —0.3 percent for p$ol Pm equal to 2.9, 
and —0.6 percent for p^ 0 lp m equal to 4.5. Al- 
though the contribution of the cold gas was itself 
influenced by the finite expansion (up to a 50- 
percent decrease at the highest pressure ratio), 
the cold gas contributes so little to laminar-flow 
attenuation at. these low pressure ratios that the 
error in assuming a negative shock is negligible. 


For higher shock pressure ratios the effect of 
the finite expansion cannot be computed. How- 
ever, examination of figure 5 indicates that, at 
a shock pressure ratio approximately equal to G, 
the contribution of the cold gas (with the negative 
shock) vanishes. At higher pressure ratios the 
cold gas tends to decrease attenuation, but this 
effect remains small compared with the contribu- 
tion of the hot gas. At a shock pressure ratio of 
20, for example, the effect of the cold-gas region 
has reached only 4 percent of the total. Thus, 
the finite expansion fan can influence only a small 
part of the total attenuation, and the assumption 
of a negative shock should give reasonably 
accurate results. 

The laminar boundary-layer finit e-expansion-fan 
solutions of reference 14 were also used to estimate 
the effect of the negative shock assumption upon 
the attenuation predict ed by reference 2 for the 
same three pressure ratios. The errors in attenu- 
ation which arise through the use of a negative 
shock arc 1.4 percent at ppjpa, equal to 1.6, 5.1 
percent at pp 0 !p m of 2.9, and 4.2 percent at 
PtiolPeo of 4.5. In this case, the contribution of 
the cold gas tends to decrease the attenuation at 
all pressure ratios of in t crest, the contribution for 
the negative shock assumption reaching a maxi- 
mum at pfoipco equal to 4.5. 

No moans of computing turbulent boundary- 
layer flow inside a finite expansion are available. 
However, according to the linearized negative 
shock approach of this paper, the eold-gas con- 
tribution is always less than 10 percent for pressure 
ratios up to 20 and always less than 10 percent for 
pressure ratios up to 10. (See fig. 5.) If the 
effect of the finite expansion in turbulent flow is 
in the same direction as in laminar flow, (lie 
cold-gas contribution, when a finite expansion fan 
is considered, would be less than that for a 
negative shock. Thus it is assumed that the 
effects of the difference between the finite-expan- 
sion and negative-shock solutions for turbulent 
flow may be neglected, at least up to shock pressure 
ratios of 10. 

The influence of nonlinear! ties in the reflection 
and transmission of waves at the entropy discon- 
tinuity was calculated for shock pressure ratios 
ppulPn equal to 1.25, 2.0, 4.0, 6.0, R.0, and 10.0. 
The nonlinearities involved required that the 
shock-tube Reynolds number and station be 
specified in order to calculate the attenuation 
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for this case. The case of aJDjv^ equal to 0.1 X 10 6 
was chosen, and the attenuations at IjD equal to 
70 and 140 were calculated for completely lam- 
inar (« = 1) flow and completely turbulent (n= 7) 
flow. 

For the laminar case, a difference of about 
one-half percent was found between linear and 
nonlinear calculations for the wave interactions 
at the entropy discontinuity for ppojp „ of 10 and 
at IjD equal to 140, the nonlinear calculation 
predicting greater attenuation. Differences at IjD 
of 70 and at lower pressure ratios were smaller and 
generally in the same direction, except at a value 
of p&ojp <» of 1.25 where the nonlinear entropy 
result gave slightly less attenuation than the 
linear result. The deviations, which were small 
in all cases, were also somewhat erratic in their 
behavior. 

Differences in the turbulent case were somewhat 
larger and ranged up to 10 percent for a value of 
IjD of 140 and a value of Pflo/p*, of 10. This 
condition is attributed to the larger relative con- 
tribution of the cold -gas region a to the total 
attenuation for turbulent flow. (See table I and 
fig. 5.) Again, the nonlinear entropy calculation 
predicted greater attenuation than the linearized 
except at a value of p&o/Pn of 1.25 where the 
trend was reversed. However, at the low pres- 
sure ratios the deviation was less than 1 percent 
and this trend is not considered to be significant. 

The effect of reflection of waves at the shock 
was estimated by the following treatment : The 
limit tenua ted shock was assumed to be overtaken 
by a single isentropic wave of strength equal to 
the total attenuation for a given condition and 
station. The resulting one-dimensional interac- 
tion was computed and the strength of the shock 
after interaction compared with the strength pre- 
dicted by the linearized theory. These results 
were computed for shock pressure ratios of 10 
and 4. 

For paoIPco equal to 10, the three cases were 
compared and the results are given in the fol- 
lowing table: 


(Too/) 

l 

Attenuation 

Percent 
difference in 

"oo 

I) 

No reflection 

With reflection 

attenuation 

0.1X10* 

140 

Li,i = — 0.60 

7i., = -0.55 

8 

.1 

140 

£*,7 = -3.57 

7-, 7= —3.16 

11 

.005 

154 

7^:= -7.0 

fi.14 

12 


The conditions of 0.005X10 fi and ^=154 

r oo H 

represent the conditions under which maximum 
attenuation would be expected for the experi- 
mental range of this paper. 

In all cases, the shock is slightly strengthened by 
the interaction, that is, attenuation is decreased. 
This is a trend which is opposite to that generally 
computed for the exact entropy discontinuity. 
For pdu/Pco equal to 4, only one case was com- 
puted, that corresponding to the largest predicted 
attenuation, a^Djv equal to 0.02X10 6 , IjD 
equal to 154. In this case, the predicted attenua- 
tion was — — equal to —2. The considera- 
P » 

tion of reflection gave — - — equal to —1.8. 

P. 

The difference in attenuation was 10 percent, 
again in the direction of decreasing attenuation. 

Tn general, then, the effects of the entropy 
discontinuity and the reflected wave at the shock 
which were neglected in the linearized theory 
tend to oppose one another and are of about the 
same order of magnitude for the worst cases 
(turbulent flow, strong shocks). 

For the pressure range of this report it has been 
shown in the preceding discussion that the two 
nonlinear effects having the largest magnitude are 
of opposite sign. The effect of neglecting the 
finite expansion-fan width is small for laminar 
flow and is also presumed to he small for turbulent 
flow. Consequently, the neglect of those three 
effects, which have a net influence much smaller 

A 

than the nonlinear P effect considered, appears 
to be justified for attenuation. Although the 
aforementioned effects arc negligible for attenua- 
tion, they can appreciably influence the pressure, 
density, and velocity distribution at points re- 
moved from the shock. This behavior arises from 

A 

the fact that the relative influence of the P waves 
of region £ decreases as the distance behind t lie 
shock increases; the influence reduces to zero at 
and behind the entropy discontinuity. (See 
pages 33 and 34 and fig. 13 of ref. 1 for further 
discussion.) 

CONCLUDING REMARKS 

The theory of NAOA Technical Note 3375, in 
which shock-wave attenuation is calculated by 
use of a linearized form of the hyperbolic system 
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of equations of motion and energy and through 
the assumption of equivalent incompressible 
steady-flow friction and heat-transfer coefficients, 
has been modified in the following manner: 

(1) Incompressible unsteady skin-friction co- 
efficients have been determined by an integral 
method. The resulting unsteady incompressible 
skin-friction coefficient is corrected for compressi- 
bility by a reference temperature method. 

(2) A nonlinear regional approach has been 
employed to permit the extensions of the theory 
to large attenuations. This approach modifies 
only the forward running waves generated in the 
hot-gas region; these waves are shown to dominate 
the attenuation process. 

(3) Transition effects are considered. The 
method assumes instantaneous transition from 
laminar to turbulent flow. The Reynolds number 
of transition then becomes a parameter of the 
attenuation problem. 

The modified theory has been evaluated for a 
range of shock pressure ratios of general experi- 
mental interest. Curves arc presented to permit 
easy prediction of attenuation for shock pressure 
ratios to 20 and a range of shock-tube Reynolds 
numbers. Results for the linearized theory with 
all laminar and all turbulent flow are compared 
with results of NACA Technical Note 3278, and 


the predicted attenuations are found to be in fair 
agreement. 

Experimental local shock pressure ratios have 
been determined for ideal shock pressure ratios 
from 1.3 to 10 for a range of shock-tube Reynolds 
numbers and position from the diaphragm station. 

Comparison of the modified nonlinear theory 
using a transit ion Reynolds number based on flow 
length of 2.5X10 6 with the experimental results 
shows good agreement except for the following 
situations: 

(1) The highest shock pressure ratios at the 
lowest shock-tube Reynolds numbers where the 
effects of the nonperfect diaphragm burst are 
believed to dominate the flow. 

(2) The lowest shock pressure ratios, for which a 
lower transition Reynolds number (1.25X10 6 ) ap- 
pears to give better agreement. 

The effects of considering a finite expansion fan, 
the exact interaction at the entropy discontinuity, 
and reflection at the shock wave (all of which are 
neglected in the present treatment) arc calculated 
for certain cases. The net effect of these three 
contributions is shown to be small compared with 
the nonlinearities accounted for in the stepwise 
regional calculation. 

Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., June 25, 1958. 


APPENDIX A 


SUMMARY OF LINEARIZED ATTENUATION EQUATIONS 

Tlie changes in the characteristic parameters P, Q, and S/R as a function of local skin friction are 
derived in reference 1 and those appropriate to shock-wave attenuation are summarized in this appendix. 
Subscripts refer to positions on the wave diagram of figure 1. 


P c -P t P-P ao 
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APPENDIX B 

DERIVATION OF EQUATION FOR Zjv 


The derivation for Z\ will be obtained and then 
generalized to arbitrary N. Tlio wave diagram of 
figure 9(a) shows the location of the numbered 
points. The symbol denotes the distance of the 
point 0 from the origin, and so forth. The follow- 
ing equations arc obtained by simple geometry 
from figure 9(a) : 

(Bl) 

? 2 = Usiti + Pn(f 7 -«i) + (P+a)i& -U) (B2) 


Combining equations (Bl) and (B2) yields 


tr. , u, h u,T-Un 

1 1 _ L ii — (Jj + a) i J 1 1 f n (P+ a ) i 


Also from figure 9(a), 

= t sib i 


(S3) 


(B4) 


T Un(t?—ti) — £ si (f 7—t\n) 

+ (P+a)i(*n-*io) (B5) 

j*u~ f'Viod - (f 7 -(“fl)i(fii fio) (Bb) 


Combining equation (B4) with equation (B6) and 
combining equation (B5) with equation (B6) 
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yields 

1 1 1 &l 

1 10 (t +u)i Uni 

(B7) 

and 

^io £ a i Un / 
U^Um-Uj \U ) 

(B8) 

Now 

/ 1 1 ^11 ^11 t IQ 

77“77 - fio^7 

(B9) 

Substitution 

of equations (B3), (B7), 

and (B8) 

into equation (B9) produces 


In 

L si-Uu a T / 

l) (BIO) 

lx 

U, t r, (C7+a),-C n \*i 

or 

Tr =z, (s _1 ) 

(Bll) 


Since the origin in figure 9 could be shifted to 
any arbitrary location and the same geometrical 
relations above could be derived in the new 
location the generalized formula for Z v can be 
written from inspection of equations (BIO) and 


(Bl 1) as 


Z, v = 


U * n — £ ;l a t +i fly 

Tj s n U n {V -f~ d ) n t 


.v+i 


(Bl 2) 
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TABLE I— ATTENUATION FUNCTIONS AND COMPRESSIBILITY CORRECTIONS 




( Gp )«, 7 

( IW..U 


fia.7 

Ha ,13 

( Hff ) ,1 

(Of)fl. 7 

( 8 p ) j3,13 

Hp.i 



(£&)«.! 

(QsO«.r 

-KHffl/s.u 

1.00 

0 

0 

0 

1. 000 

1.000 

1.000 

0 

0 

0 

1.000 

1.000 

1.000 

0 

0 

0 

1.25 

-.0367 

-.000223 

-.0000410 

.996 

.980 

.976 

- 2. 627 

-.0142 

-.00255 

1.004 

1.018 

1.022 

- 2 . 663 

-.0144 

-.00259 

1.50 

-.0844 

-.000943 

-.000200 

.992 

.962 

.954 

- 3. 775 

-.0338 

-.00695 

1.006 

1.032 

1.038 

- 3.859 

-.0348 

-.00715 

1.75 

-.126 

-.00200 

-.000401 

.089 

.945 

.933 

- 4.705 

-.0549 

-.0121 

1.009 

1.043 

1.052 

- 4 . 832 

-.0569 

-.0126 

2.0 

-.161 

-.00326 

-.000795 

.986 

.928 

,914 

- 5. 531 

-.0765 

-.0176 

1.011 

1.052 

1.064 

- 5.692 

-.0797 

-.0184 

2.4 

-.199 

-.00531 

-.00138 

.981 

.903 

.884 

- 6. 733 

-.111 

-.0269 

1.014 

1.065 

1.080 

- 6 . 932 

-.116 

-.0282 

2.6 

-.210 

-.00626 

-.00167 

.978 

.891 

.870 

- 7. 302 

-.128 

-.0316 

1. 015 

1.071 

1.087 

- 7. 512 

-.135 

-.0333 

3.0 

-.219 

-.00785 

-.00218 

.974 

.868 

.843 

- 8.394 

163 

-.0412 

1.017 

1.081 

1.099 

-8 613 

-.171 

-.0434 

3.6 

-.205 

-.00919 ] 

-. 00267 

.967 

.835 

.804 

- 9.961 

-.214 

-.0558 

1.020 

1.094 

1.115 

- 10. 166 

-.224 

-.0584 

4.0 

-.180 

-.00925 

-.00275 

.963 

.815 

.780 

- 10.973 

-.248 

-.0655 

1.022 

1.102 

1.125 

- 11.153 

• 257 

-.0683 

5.0 

-.081 

-.00637 

-.00199 

.953 

. 165 

.723 

- 13.444 

-.332 

-.0900 

1.025 

1.118 

1.145 

- 13. 525 

-.339 

0920 

G.O 

.056 

.GOODS 

.00003 

.943 

.719 

.671 

- 15. 860 

-.416 

-.114 

1.028 

1. 131 

1.161 

- 15.804 

-.415 

-.114 

( 

7.0 

.214 

.0115 

.00370 

.933 

.675 

.622 

- 18.243 

-.498 

-.136 

1.031 

1.142 

1.175 

- 18.029 

-.486 

-.132 

8.0 

.386 

.0255 

.00863 

.924 

. 634 | 

.576 

- 20.603 

-.580 

-.163 

1.033 

1.152 

1.187 

- 20.217 

-. 555 

-.154 

9.0 

.560 

.0418 

.0144 

.914 

. 596 ; 

.534 

- 22 . 953 1 

-.661 

-.187 

1.035 

1.160 

1.197 

- 22.393 

-.020 

-.173 

10.0 

.736 

.0600 

.0210 

.905 

.559 

.494 

- 25 . 292 

-.743 

-.211 

1. 036 

1. 168 

1.207 

- 24.556 

-.683 

; -.190 

12.0 

1.072 

.0995 

.0357 

.886 

.490 1 

.422 

- 29.945 

-.904 

-.259 

1.039 

1.180 

1.222 

- 28.873 

-.804 

-.224 

15.0 

1.514 

.160 

.0592 

.859 

.401 1 

.330 

- 36.911 

- 1.143 

-.331 

1.043 

1.194 

1.240 

- 35. 397 

-.983 

-.271 

20.0 

1.756 

.213 

.0805 

.812 

.278 

.212 

- 49. 167 

- 1.559 

-.454 

1.047 

1.211 

1 . 2 G 1 

- 47. 411 

- 1.346 

1 -.373 
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